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Abstract 

The AMS-02 has measured the cosmic ray electron (plus positron) spectrum up to ~TeV with 
an unprecedent precision. The spectrum can be well described by a power law without any obvious 
features above 10 GeV. The satellite instrument Dark Matter Particle Explorer (DAMPE), which 
was launched a year ago, will measure the electron spectrum up to 10 TeV with a high energy 
resolution. The cosmic electrons beyond TeV may be attributed to few local cosmic ray sources, 
such as supernova remnants. ‘Therefore, spectral features, such as cutoff and bumps, can be 
expected at high energies. In this work we give a careful study on the perspective of the electron 
spectrum beyond TeV. We first examine our astrophysical source models on the latest leptonic data 
of AMS-02 to give a self-consistent picture. Then we focus on the discussion about the candidate 
sources which could be electron contributors above TeV. Depending on the properties of the local 
sources (especially on the nature of Vela), DAMPE may detect interesting features in the electron 


spectrum above TeV in the future. 
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I. INTRODUCTION 


The Alpha Magnetic Spectrometer (AMS-02) launched in May 2011 has taken the mea- 
surement of cosmic ray (CR) leptonic spectra to a new level [1]. Unprecedentedly precise 
results have been released in the energy range from ~ 1 GeV to ~ 500 GeV for intensities of 
both electron and positron. The well-known electron/positron excess, which was uncovered 
by earlier satellite experiments such as the Payload for Antimatter Matter Exploration and 
Light-Nuclei Astrophysics (PAMELA) Ld l and the Fermi Large Area Telescope (Fermi- 
LAT) |4, |5], and also the anterior balloon-borne Advanced Thin Ionization Calorimeter 
(ATIC) experiment |6], has been confirmed by AMS-02. Sources which have the potential to 
provide primary electron/positron pairs, e.g., pulsars and annihilation of dark matter (DM), 
are thus involved into CR models to interpret those excesses. However, some global fittings 


of AMS-02 leptonic spectra indicate that the electron spectrum has a larger excess to the 


background than that of positron . As the contributions of e* pairs from exotic sources 
like pulsars or DM are constrained by the positron spectrum, these sources seem not enough 


to explain the total electron excess. 


A possible interpretation of the extra electron excess is the hardening in high energy 
range of the electron spectrum of the supernova remnant (SNR) background, which can be 
attributed to the fluctuation given by local discrete SNRs |9]. Thus the concept of dividing 
local SNRs and distant SNRs steps again into the spotlight. This concept was first put 
forward by Shen [10] and improved by later works. In the model of Atoyan et al. m 
a nearby (< 100 pc) and relatively young (< 10° yr) source and continuously distributed 
distant sources (> 1 kpc) contribute separately to the electron spectrum; they also adopted 
a energy dependent diffusion coefficient in the propagation model. Kobayashi et al. 
went further on this scenario by using several real sources with known ages and distances 
as local electron accelerators. After the publishing of AMS-02 data, Di Mauro et al. |13 
fit all the leptonic data simultaneously, applying the method similar to the works above 
when dealing with local and distant SNRs. They derived spectral index and normalization 
of injection spectrum of individual SNR from radio observations, comparing with uniform 


injection spectral parameters of all the sources adopted in Ref. id 


In fact, there are alternative explanations toward the electron/positron excess, such as 


improper propagation parameters used in previous works |14. In these cases, exceptional 
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consideration of local sources may not be necessary to keep consistency with the AMS-02 
data. However, F Moa based Cherenkov telescopes like the High Energy Stereoscopic 
System (HE 13 id and the Very Energetic Radiation Imaging Telescope Array System 
(VERITAS) ph 1T seemed to detect a cut-off around ~ 1 TeV in the e~ +e" spectrum, which 
cannot be described by a continuous distributed background. The Dark Matter Particle 
Explorer (DAMPE) P 


of 5 GeV—10 TeV with unprecedented energy resolution (1.5% at 100 GeV). As nearby 


launched in December 2015 aims to measure electrons in the range 


sources have the potential to contribute to the highest energy range covered by DAMPE, we 
can expect to see some spectral features in future DAMPE results. In this case, separating 
local sources from the continuous distribution would be inevitable. 

Basing on previous works mentioned above, we perform a careful analysis of the local 
SNRs and their parameters of injection spectra. We assume that pulsars are extra positron 
sources, and perform global fittings to the latest leptonic data of AMS-02 for several SNR 
parameter settings. We show below that the choice of parameters of a particular SNR 
has a significant influence on its contribution. Although the electron energy range covered 
by AMS-02 is under TeV, fittings to the AMS-02 leptonic data provide a self-consistent 
picture for the astrophysical source models. As the local sources accounting for the AMS-02 
results may provide contribution to the TeV scale, the AMS-02 data could also constrain 
the properties of the predicted e~ + e* spectrum above ~ TeV. Combining with the fitting 
results, we then discuss the parameters of the local sources which have the potential to 
contribute to TeV and give further predictions of the e~ + e* spectrum up to the energy 
range of 10 TeV, which can be measured by DAMPE. 

This paper is organized as follows. In Sec. we describe our calculation towards the 
injection and propagation of Galactic electrons and positrons to get leptonic spectra. The 
results of global fittings to leptonic data of AMS-02 and our further predictions to the 
electron spectrum in the TeV range are presented in Sec. [M] and Sec. respectively. We 


summarize our work in the last section. 


Il METHOD 


In this section, we introduce the semi-analytical solution to the propagation equation in 


the first part. Then we discuss the possibly that SNRs are the most important sources of 
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high energy CR electrons. Discussions of other sources, including pulsars and secondary 


electrons/positrons, are given in the later subsections. 


A. Propagation of Cosmic Rays 


The propagation of CR electrons in the Galaxy can be described by the diffusion equation 


with additional consideration of energy loss during their journey, which may be written as 
< — V(DVN) - (bE) =, (1) 


where N is the number density of particles, D denotes the diffusion coefficient, b denotes 
the energy-loss rate and Q is the CR source function. Galactic convection and diffusive 
reacceleration are not taken into account here, since they have little effect above 10 GeV 
19]. We treat the propagation zone of CRs as a cylindrical slab, with radius of 20 kpc [19 
and a half thickness z}. The diffusion coefficient D depends on the energy of CRs, which has 
the form D(E) = BDo(R/1GV)°, where Do and 6 are both constants, 8 is the velocity of 
particles in the unit of light speed and R is the rigidity of CRs. To give a constraint of major 
propagation parameters—(Dp,0,z,), Boron-to-Carbon ratio (B/C) is widely used. Unstable- 
to-stable beryllium nae /°Be) is also helpful to constrain CR propagation. We adopt the 
B/C af of sce pi , AMS-02 aL 10Be/°Be data of Ulysses [2 mt ACE [23], Voyager l, 
IMP l, ISEE-3 25), ISOMAX [2 26), and embed the CR propagation code in the Markov 
Chain Monte Carlo (MCMC) sampler to acquire best fitted propagation parameters (this 
work is in preparation). We adopt Do = 2.12 x 1078 cm?s~!, 6 = 0.548, zn = 3.8 kpc in this 
work. 

Positrons and electrons with energy higher than 10 GeV suffer from energy loss during 
their propagation mainly by synchrotron radiation in the Galactic magnetic field and inverse 
Compton radiation in the interstellar radiation field consisting of stellar radiation, reemited 
infrared radiation from dust, and cosmic microwave background (CMB). We set the inter- 
stellar magnetic field in the Galaxy to be 1 uG to get the synchrotron term [19]. For the 
inverse Compton process, if we use the cross section for Thomson scattering, the energy-loss 


rate b(E) has a quadratic dependence on energy: 


-= (2) 
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FIG. 1. Energy-loss coefficient bọ as a function of e* energy. The dashed line stands for the 
case of Thomson approximation and the solid line describes the drops of bg due to Klein-Nishina 


correction. 


where bo is a constant. However, when the energy comes higher, a relativistic correction 
to the cross section, namely a Klein-Nishina cross section, is needed. Here we adopt the 
description of Ref. which gives a reconcilable approximation between Thomson limit 
and Klein-Nishina limit. The temperature and energy density of radiation field components 
are: 20000 K and 0.09eVcm~? for type B stars, 5000 K and 0.3eVcm~? for type G-K 
stars, 20 K and 0.4eV cm ~® for infrared dust, and 2.7 K and 0.25 eV cm~? for CMB. In this 
case, bp is not a constant anymore but decreases with the energy as shown in Fig. [| We 
can find from Fig. [I] that the relativistic correction becomes important as long as energies 
of electrons/positrons are higher than 10 GeV. We still use the symbol bọ, which has the 


connotation of bo(£), in our later work for convenience. 


Cylindrical coordinate is applied in our work to describe the disc-like geometry of the 


propagation zone, and the location of the Earth is set to be zero. For a point source with 
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burst-like injection, the source function can be written as 
Q(E,t,r, z) = Q(E)ó(t — t,)6(r — r5)6(z — 2s) , (3) 


where Q(E) represents the energy distribution of injection, ts is the time of CR injection, 
rs and z, are radial and verticle location of the source, respectively. As long as Galactic 
CR sources are mostly distributed in a much thinner vertical scale comparing with z, and 
so does the solar system, we assume z, = 0 in our work for all the sources. Then the 
time-dependent Eq. (@) can be solved semi-analytically with the help of Green’s function 


G(E,t.r + E,,ts,r;) working in the Fourier space. We follow the Green’s function used in 


Ref. id. 


G(E,t,r + E,, ts, rs) = 6(E, — Eo) = - = |- = 


: (4) 
20 Dy 
rae): 
where kn = (2n + 1)a/(2z%,), Eo = E/[1 — bo E(t — t,)], and 


is the diffusion distance for particles with initial energy Æo and final energy E. The solution 


of Eq. has the form of G(E,t,r © E,,ts,7,)Q(Es, ts, 1s), so the observed CRs contributed 


by a source with distance r and age t should be expressed as 


1 E r? E 
N(E,t,r) = TA (1 = bo Et) ? exp (-5) Q (r) 


«ieo (2H) 7 


However, considering the efficiency of doing numerical calculation, Eq. (6) may not be a good 


expression since we need to include more terms in higher energy range to guarantee precision 


of the calculation. Thus a spherically symmetric time-dependent solution to Eq. (1), which 


has a form of bs 


N(E,t,r) = reas (1 — bo Et)’ exp (-5) Q (r) , (7) 


can be treated as a substitute of Eq. (6). Given the disc-like geometry of propagation zone, 


Eq. (Z) is valid when A < zp. Assuming F is close enough to Eo, the diffusion distance 
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FIG. 2. Comparision between diffusion distance of CRs and half-thickness of propagation zone. 


The solid line represents the diffusion scale A ~ 2\/D(E)t while the dashed line marks zp. 


defined in Eq. can be approximated by 2,/D(E)t. As shown in Fig. P| the condition 
À < Zp, is always satisfied as long as Æ > 10 GeV, and the difference increases in higher 
energy. In fact, we find the error is less than 1% for Æ > 10 GeV after calculating ratios 


between these two expressions. 


B. SNRs as Electron Source 
1. Local SNRs and SNR Background 


SNRs are believed to be the main astrophysical sources of primary Galactic CRs, such as 
nucleus and electrons. Particles can be boosted to very high energy through diffusive shock 
wave acceleration in SNR. However, among the accelerated particles, electrons/positrons 
undergo significant energy loss during their propagation through electromagnetic radiation. 
Eq. indicates the life time of an electron is roughly 1/(b)£), thus electrons with higher 
energy become inactive faster. Since ô < 1, the diffusion distance 2,/D(E\t has an inverse 
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relationship with electron energy. For example, electrons in ~ TeV fade within a radius of 
roughly 1 kpc from their source. Thus high energy part of electron spectrum can only be 
contributed by several local sources, and a simple continuous source distribution is no longer 
valid due to the spectral fluctuations induced by those few sources. Then it is important 
to separate local discrete sources from distant sources in calculation, as first proposed by 
Shen [10]. We treat SNRs within 1 kpc as local sources and farther sources as background 
contributors of electrons la. The intensity of electrons from a SNR nearby can be simply 
expressed by 


I(E) = (Et 4) (8) 


In order to obtain distant components, we calculate the electron spectrum produced by 
a smooth distribution of SNRs in whole range of distance and age first, and then sub- 
tract the local components in continuous form. We set supernova explosion rate to be 
f = 4 century! galaxy! Ri Since about 2/3 of supernovae are expected to be type II 
supernova, we can use the population of pulsars to describe the spatial distribution of SNRs. 


Here we choose the Galactic distribution of pulsars given by Ref. [30] as 


p(R) ~ R” exp (-=) (9) 


where n = 2.35,0 = 1.528 kpc and R is the distance to the Galactic center. Note that the 
zero point of this distribution is the Galactic center, rather than the solar system used in 
our work. Thus p should also depend on the azimuth angle y. Here we aim to find how 
local sources create spectrum features in TeV range. As the diffusion distance of 1 kpc 
corresponds to a electron cooling time of roughly 3 x 10° yr, sources older than this age are 
treated as background SNRs in our calculation. Finally we get the electron spectrum of 


background component: 


iE foe) 2m tm Fin 2m 
Ds | a f ar | d- f a f ar f dip 
4T 0 0 0 0 0 0 (10) 
x fpr, p) N(E,t,r)r, 


where p(r, p) is the normalized distribution, rm = 1 kpc and tm = 3 x 10° yr for the case 
1/(boE) > 3 x 10° yr (otherwise tm takes 1/(bp£)). Indeed, full propagation equation with 
consideration_of convection and reacceleration may be solved with public numerical tool 


GALPROP which can give a more accurate result. Nevertheless, the spatial zero point 
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TABLE I. Basic parameters for local SNRs within 1 kpc. The columns show the Green catalog 
name, the association name, the radio flux BGHZ, the radio spectrum index a,, the distance r and 


the SNR age t. Note that t is the observed age, the actual age should be T=t+r/c. 


used in GALPROP is the center of the Galaxy, we may need to do some troublesome work 
to divide local sources from the background. That’s the reason we choose an semi-analytic 


treatment toward Eq. (i). 


2. Parameters of SNRs 


For SNRs, particles are accelerated by shock acceleration mechanism, or say, the first 
order Fermi acceleration. The energy spectrum produced by Fermi acceleration is thought 
to be a power law form. Taking energy loss and escape of particles into account, the emergent 


spectrum of SNR can be described by a power law form with an exponential cutoff: 
Q(E) = Qo(E/1 GeV) 'exp(—E/E,) , (11) 


where Qo is the normalization of the injection spectrum. Note that distant SNRs are treated 
as a continuous distribution. Thus we assume that they share common Qo and y, and set 


these two parameters to be free in following fittings. The energy cut-off is fixed at 20 TeV for 
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background sources. For local SNRs, we attempt to investigate their parameters individually. 
Table [I] lists objects locating within 1 kpc included by the Green catalog of SNRs (of 
course those with a measured distance) with two additional sources Monogem Ring and 
Loop I. 

Several SNRs have gone through multi-wavelength measurements, from radio to y-ray 
bands, which are helpful to estimate parameters in Eq. (II). These sources include HB9 
(radio and y-ray), Vela Jr. (radio, X-ray and y-ray), RX J1713.7-3496 (radio, X-ray and 7- 
ray), and Cygnus loop (radio and y-ray). Note that we just pick observations with available 
data. Radio and X-ray emissions are produced by the synchrotron radiation of relativistic 
electrons accelerated in the SNR, while y-ray emissions could have either a leptonic origin 
or a hadronic origin, which are generated through scattering of background photons by 
relativistic electrons or through 7° decay resulting from the collision of accelerated protons 
with ions in the background plasma. As we are interested in the electron spectrum, a purely 
leptonic model is the priority in our fitting. If the leptonic fitting fails, this model is replaced 
by a hybrid model, where contributions from electrons and protons are comparable in y-ray 
spectrum. Purely hadronic model would not be discussed in our work. 

For the leptonic model, the energy spectrum of accelerated electrons still has the form 
of Eq. (11), while we remark the parameters in those expression as Qoe, Ye and Ece, to 
distinguish from the case of proton. The background radiation fields consist of interstellar 
infrared radiation, optical radiation and CMB. The contributions to the inverse Compton 
y-ray spectrum from the interstellar radiation field (ISRF) are more important than that 
from CMB . Here we adopt the ISRF model given by Ref. . In the hybrid model, the 
energy spectrum of protons has the same form as that of electrons, where Qop, Yp and Eep 
are corresponding normalization parameter, spectral index and cut-off energy of protons. 
Assuming charged particles share the same acceleration mechanism, the spectral index of 
protons could be identical to that of electrons lea eh, Thus the four free parameters in 
leptonic model are Qo, Ye, Ece and magnetic field B. For the case of hybrid model, there 
are three additional parameters, Qop, Ec,p and number density of target proton Ny. 

We apply naima, a python-based package, for computing non-thermal radiation processes 
and fitting the spectral energy distributions kd. It uses Markov-Chain Monte Carlo emcee 
samptine od to find best-fit parameters of physical radiation models and thus determines 


the radiation mechanism behind the observed emission. The parameters of the best-fit 
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FIG. 3. Energy spectrum distribution of the sources Cygnus Loop, HB9, Vela Jr. and RX J1713.7- 
3946. The components are: synchrotron (green line), IC(blue line), non-thermal bremsstrahlung 
(purple line), neutral pion decay (red line). Total y-ray emission for hybrid scenario is marked by 


yellow line. References for the fitted data points are given in Table [I] 


models are compiled in Table [Hand the best-fit spectra are shown in Fig. 

Different from four SNRs mentioned above, we need other approaches to determine pa- 
rameters for other SNRs without precise high energy y-ray observations. Assuming the radio 
emissions of SNRs are entirely induced by synchrotron emissions of electrons, Di Mauro et al. 

13] provided an estimation method of Qo as 


m F es B —(y+1)/2 BiGHe 
= 12x 104 Gev— (0.79)7{ —— | {| —2— r 12 
AE (i) (samc) ( J ). (12) 


which relies on spectral index y, source distance r, magnetic field B in SNR and the radio 
flux at 1 GHz B19". For the synchrotron emission of electron system, the electron spectral 
index can be derived from the photon spectral index. Thus we have y = 2a, + 1, where a, 
is the radio spectral index. We list B!C#* and a, in Table [] Another critical parameter B 
can only be approximately estimated by several methods, such as Zeeman effect, Faraday 
rotation, and minimum energy or equipartition. The equipartition or minimum energy 
method calculates magnetic field only depending on the radio synchrotron emission of source 


. It is a useful tool if no other data of the source is available. According to the result of 
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Arbutina et al. lg, we can obtain the magnetic field of SNR as 


1 (TEAT 
BIG] = [6.286 x 1007779/2 ae Mak a a) i ) me 
aa MCP) (13) 
2c, )A-1/2 BiGHzy 2/(y+5) 
ye) (1+4) T. y] 
C5 f r{kpe] 6[arcmin]? 


where mec? is rest energy of electron, cı and cs = csl' (24+) (5*)/(7 + 1) are defined in 
Ref. [70], f is the volume filling factor of radio emission, 0 is the angular radius of SNR 
which can be found in Green catalog. In Eq. (13), « is a parameter depending on y and 
ion abundances of SNR. We adopt a typical value f = 0.25 and simple ISM abundance 
H:He=10:1 in our estimation. Note that this minimum energy method is applicable only 
for mature SNRs with a > 0.5. For G114.3+00.3 and G127.1+00.5, we take a typical value 
of B = 30 uG. Finally, We follow the method proposed by Yamazaki et al. mi to give an 
estimation to Ee. When the age of a SNR meets t > 10° yr, synchrotron loss restricts the 


maximum energy of electrons. Thus we have 


E. = 14 TeV ho? v, 8 (B/10uG)-"?, (14) 


1 which depends on evolution of SNR. 


where vs s is the shock velocity in unit of 108 cm s7 
See details in Ref. for the definition of h and here we take it to be unit. 

So far we have estimated parameters for all the sources in Table [except for Vela(XYZ), 
Monogem Ring and Loop I. The results are listed in Table] Vela(XYZ) is generally believed 
to be an important local source of CRs. We discuss it in detail in the next section. For 
Monogem ring and Loop I, they are classified as possible or probable SNRs by Green and 
are not included in the catalog due to the lack of understanding to them. We treat them 
as potential sources of electrons in the next section and set their parameters to be free. 
Note that all those estimations of SNR parameters rely on the electromagnetic radiations 
of SNR. These emissions are not produced by electrons observed today but some ’younger’ 
ones in SNRs. This implies that even if observations of photon emission are precise enough, 
the derived electron parameters may be different from those of injected electrons due to 
the variation of parameters along with the evolution of SNRs. Therefore the aim of our 
estimation is to give some reference values for electron injection parameters, but far from to 


‘determine’ them. 
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TABLE II. The electron spectrum parameters for our sample. The parameters of HB9, Vela Jr. and 


RXJ1713.7-3946 are obtained by fitting the multi-wavelength emission by leptonic model. The hy- 
brid model give better fitting for Cygnus loop, and the extra parameter are Qo,» = 1.5x 10°? Gey, 
Ecp = 12GeV, Ny = 0.46cm~3. Magnetic field strength of G065.3+05.7 and G156.2+05.7 are 
deduced from Eq. (13), and we adopt 30 uG for G114.3+00.3 and G127.1+00.5. For those sources 
Qo,e and Eee are derived from Eq. and Eq. respectively, while ye is given by the relation 


y = 2a, +1. References for the multi-wavelength fitting data are given in the last column. 


C. PWN as Electron and Positron Source 


Pulsars are known to be the most_important astrophysical sources of high energy elec- 
tron/positron pairs in the Galaxy |85]. They produce relativistic wind carrying charged 
particles at the cost of their spin-down energy kd. Since pulsar is formed in SN explosion, 
it is initially surrounded by its companion SNR. When the relativistic pulsar wind impacts 
on the cold SN ejecta which expands with a slower velocity, a termination shock is formed 
besides a forward shock. The termination shock propagates inward and reaches the radius 
where the outward pressure of pulsar wind balances the internal pressure of the shocked 
bubble (see Ref. 


wind nebula (PWN). After particles inter the PWN, they are confined by the magnetic field 


and references therein). The shocked region is dubbed the pulsar 


here for a long period, until the crushing of the PWN. Thus the spectrum of electrons and 
positrons injected into ISM should be the spectrum inside the PWN when it is disrupted, 
other than the spectrum inside the magnetosphere of pulsar ld, Like SNRs, pulsars can be 
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divided into local pulsars and smooth distant components. Delahaye et al. find that the 
contribution from pulsar background is negligible compared to those from local ones, thus 
we do not take the former into account in our calculation. Parameters of nearby pulsars can 
be found in the ATNF catalog ; 

We assume the injection e* spectra of PWNs have the same form of Eq. (I); a burst-like 


injection spectrum is also adopted. Note that the spectral index used here is associated with 
the PWN, and cannot be derived from the spectral index of the pulsed radio emission from 
pulsar given in the ATNF catalog. Thus if the radio spectral index of PWN is not available, 
we set Ypwn as a free parameter in the following fittings. The normalization parameter Qo 


is linked to the spin-down energy W, dissipated by pulsar by: 


f QEEaE = nw, (15) 


where Emin = 0.1 GeV and ņ is the efficiency of energy conversion treated as another free 
parameter. The spin-down luminosity E evolves with the age of pulsar t as E = E,(1 + 
t/T)~? 89 , where Æp is the initial spin-down luminosity, 7 is the spin-down time scale of 
the pulsar, assumed to be 10 kyr in our work. Integrating E with time, we get the expression 
of spin-down energy Wp = Et (1+t/7), where E and t can be found in ATNF catalog. For 
the cut-off energy, We take Ee = 2 TeV following the work of Ref. . 


D. Secondary Electrons and Positrons 


Secondary electrons and positrons are created by inelastic collision between CR. nuclei 
and ISM. The CR nuclei mainly consist of protons and a particles while H and He are major 
components of ISM. As can be seen from the AMS-02 result, the positron fraction is smaller 
than 10% below ~ 100 GeV, where secondary positrons should contribute less than 10%. 
Since the spallation process produces more positrons than electrons, secondary electrons 
possess a further smaller percentage comparing to the total electron intensity in the range 
mentioned above. Thus we neglect the secondary electron component and concentrate on 
secondary positrons. In this part, our calculation follows the method of Ref. 00 . The source 


function of secondary e* is assumed to be steady and homogeneous in a slab geometry: 


dol BLE 
Ga a f dE! ®;(E") fol) (16) 
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so that the propagation equation can be solved semi-analytically with a relatively simple 
form (see for details). In Eq. (46), ¿ and j mark the species of CR nuclei and ISM gas 


respectively. The number density of ISM is set to be ny = 0.9 cm~? and nye = 0.1 cm~’. The 


. Di Mauro 


intensities of incident CR nuclei are donated by ® which can be estimated by the observed 
intensities at earth. We employ the expression of ® proposed by Shikaze et al. 

et al. fit the AMS-02 data of proton and Helium to refresh the parameters in the model 
of ®. For differential scattering cross-section do/dE, Kamae et al. 02 provides functional 
formulae for proton-proton collision. Empirical rescaling based on this result are applied 
to estimate cross-section of collision between other species 03). Besides, we introduce a 
free parameter Ce+ to rescale the secondary intensity, considering the uncertainty in the 


calculation above, to accommodate the data. 


Ill. FITTING TO THE AMS-02 DATA 


Up to now, AMS-02 provides the most precise measurement of leptonic spectra below 
1 TeV. We attempt to perform global fittings to all the leptonic data released by the AMS- 
02 Collaboration, including the positron fraction, positron plus electron spectrum, positron 
spectrum and electron spectrum [94H96]. The global fitting is useful to set stringent con- 
straints for astrophysical contributors, and leads to a self-consistent picture 

Astrophysical components considered in our model have been discussed in the previous 
section, including background SNRs as main contributors of electrons in low energy range, 
local SNRs as dominant contributors of electrons in higher energy range, secondary positrons 
which dominate low energy part of positron spectrum, and local PWNs as predominant 
components in higher energy range of positrons. 

In this section, we explain the AMS-02 data by using several astrophysical source models. 
The characteristic of each model depends on which sources are chosen to be predominant 
local SNRs. We can see below that according to the calculation in the previous section, local 
SNRs listed in Table] hardly have significant contributions to electron intensity within the 
energy range of AMS-02. They would not play the role of predominant local sources in this 
section; we simply add their contributions for each model. 

Positrons contribute only ~ 10% of the total e~ + e* intensity. We would like to simplify 
the constitution of positron spectrum, and use a single PWN to fit the high energy part of 
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e*t spectrum. Di Mauro et al. have already given a ’single-source’ analysis in their work 
and find that Geminga is the most proper one among their candidates. However, due to the 
old age of Geminga (342 kyr), a spectral roll-off may appear below 1 TeV. Thus if Geminga 
is chosen as the single PWN, it may induce a slight break in e~ +e* spectrum below 1 TeV. 

Another famous pulsar B0656+14, also namely Monogem, is believed to be associated 
with Monogem Ring. It lies at a distance of 0.28 kpc with an age of 112 kyr which is younger 
than Geminga |88]. Monogem has the potential to significantly contribute to the positron 
spectrum from 100 GeV to 1 TeV (see Table 4 of Ref. ). Also, we can see from Figure 2 
of Ref. Ri the spectrum of Monogem has a similar shape with that of total PWNs in the 
ATNF catalog. Therefore we expect a PWN with the similar distance and age as Monogem 
can play a role as the single positron source. 

We would show in following subsections that Monogem does well in the fittings as a single 
positron source. Comparing with the total spin-down energy of Geminga 1.26 x 10*° erg and 
the required ņ of only 0.27 in the ’single source’ model in Ref. ha, the fittings for Monogem 
requires a large 7 (0.6~0.8 in our models) due to its low spin-down energy of 1.78 x 10* erg. 
However, Monogem can extent its electron/positron spectra to higher energies than Geminga 
because of its smaller lifetime, and is helpful to explain the cut-off indicated by HESS and 
VERITAS data. In this work, we take Monogem as the single positron source in our fittings. 

Hence, for different models shown below, the common free parameters are Qo and y for 
SNR background, rescaling parameter c.+ for secondary positron, 7 and ypwn for Monogem 


and a solar modulation potential ¢ to accommodate data below tens of GeV. 


A. Vela YZ 


It is widely believed that the famous and well studied SNR Vela is an important local 
source of Galactic CRs due to its appropriate age and distance and its strong radio emission 
12), |13]. Fig. [4]shows the contours of electron intensity as a function of r and t of source at 
different energy, assuming typical input energy (10*° erg), spectral index (2.0) and cut-off 
energy (10 TeV) of electrons. Every source in Table [I| is marked in Fig. Æ| It is clearly 
that Vela predominates over other local SNRs above hundreds of GeV, and this is where the 
observed electron excess to continuous SNR model appears. 


Vela shows a shell-like radio structure consisting of three principle regions, dubbed Vela X, 
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FIG. 4. Contours of electron intensity as a function of ages and distances of local SNRs. Top 


left: Æ = 100 GeV; top right: Æ = 500 GeV; bottom left: E = 1TeV; bottom right: Æ = 5 TeV. 


Contour lines in each graph represent for E°7(E£) = 10~', 1072, 1073, 1074 GeV? cm~? s7} sr™t}, as 
noted in each graph. In the calculation, the spectral index, cut-off energy and total input energy 


048 erg, respectively. Sources listed in Table [I] are 


of electrons are assumed to be 2.0, 10 TeV and 1 
also marked in each graph. Since every SNR is constrained to share the same injection spectrum 


here, this figure is for reference only. 


Y and Z los. A weaker component Vela W observed by Alvarez et al. is not considered 
in our work. Milne 09 find the spectrum of Vela X was remarkably flat than that of Vela Y 
and Z. This unusual spectrum was first explained by Weiler and Panagia ho. They argued 
that Vela X should belong to plerions and PWN s like Crab, rather than shell-type SNRs like 
its siblings Vela Y and Z. This point of view has been accepted in later works. Consequently, 
when we estimate the contribution from Vela SNR, we ought to divide Vela X from Vela YZ. 
PWNs have different mechanisms of acceleration and evolution from shell-type SNRs. This 
means that they cannot share parameters, such as spectral index or cut-off energy, with the 


later. 


In the first model, we set Vela YZ as the unique predominant SNR and aim to check 
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if it can give a good interpretation to AMS-02 data. As suggested by Ref. , we fix the 
magnetic field to be 30 wG and cut-off energy to be 2 TeV. In the Green catalog of SNRs, the 
radio spectral index of Vela is denoted as ’varies’ perhaps due to the discrepancy between 
Vela X and YZ. We put spectral index Quela of Vela YZ into the group of free parameters 
in this model. A free Aye1q leads to the uncertainty of B!S"%. The radio flux of Vela YZ at 
960 MHz is measured to be 1100 Jy sa) so we set B1G4# = 1000 Jy here as an estimation. 
We seek the best-fit model by minimizing chi-squared statics between model and data points. 
The result of global fitting to AMS-02 data are shown in Fig. [5] In each sub-graph, all the 
components are drawn to show their contribution and black solid line represents the fitting 
result. The best-fit reduced x? is 0.473 for 182 degrees of freedom (d.o.f.), while the best-fit 
parameters are compiled in Table Solar modulation @ converges to zero in this case. 
As can be seen that this model well explains the AMS-02 data. The best-fit radio spectral 
index of Vela YZ seems reasonable, as it is close to the typical value 0.5. Di Mauro et al. [13 
also fit the four leptonic observables of AMS-02 simultaneously and obtain fairly good result. 
However, our model has two main differences from theirs. First, the propagation parameters 
they adopted is the MED model proposed by Ref. which is based on the B/C analysis 
performed in Ref. , and the nuclei data behind is taken from earlier balloon and space 
experiments; as described in Sec. we include the latest B/C data from AMS-02 in our 
calculation of propagation parameters. Besides, Vela X and YZ are assumed to provide a 
whole SNR contribution in their fittings. As mentioned above, these objects may not share 
the same spectral index, cut-off energy, and normalization parameter Qo. If the contribution 
of Vela X is considered in the energy range covered by AMS-02, it should have contributed 
to positron intensity due to its PWN nature, and would not induce a spectral structure at 


high energies above TeV. 


B. Vela YZ + Monogem Ring 


In this model, we attempt to investigate the parameters of Vela YZ at first, rather than 
taking typical values as in the previous model. The geometry of Vela SNR can be sketched 
by two hemisphere with different scales due to the asymmetrical density distribution of its 
surrounding medium [|103]. Vela Y and Z are located in the north-east(NE) part. “on had 


the equilibrium between thermal pressure and magnetic pressure, Sushch and Hnatyk [103 
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TABLE III. Fitting results of models presented in Sec. Loop I is abbreviated to Lp1 in this table. 
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FIG. 5. The results of global fitting to AMS-02 data (Vela YZ model). Top left: electron intensity; 
top right: positron intensity; bottom left: positron plus electron intensity; bottom right: positron 
fraction. The legend explains different components in each sub-graph. In the legends, TOT’ stands 
for total value, ’BKG’ stands for SNR background, ’Other SNRs’ refers to the summation of SNRs 
in Table II. Besides AMS-02, data from HESS and VERITAS measurements are also shown in 


et +e7 graph. 


give an estimation of Byp = 46 uG with the formula Byg = /8rnnekpl ye, where density 


nyng and temperature Typ are derived by Ref. . We still cite Eq. (4) to calculate cut- 


1 1 


off energy. But instead of taking typical value 10°cms~', we choose vyg = 6 x 10’cms— 
estimated by Ref. [104]. The principle behind this estimation is the expression of shock 
radius as a function of the age of SNR given by Ref. [105]. We get Ee = 4 TeV. For a, of 
Vela YZ, Dwarakanath claims a value of 0.53 combining his 34.5 MHz observation and 
other observations up to 2700 MHz. The latest work about radio spectrum of Vela is done 
in Ref. [51]. The authors of Ref. included more data points at other frequencies while 
subtracted the flux density at 34.5 MHz for the probable absorption of this measurement. 


The a, of Vela YZ is derived to be 0.735 in this work. For Qyejq = 0.53, Boo is derived to 
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FIG. 6. Electron spectra of Vela YZ with different parameters. The solid line describes the case 
B= 30 uG, Ee = 2 TeV and Quela = 0.519. The first two are typical values and the last is given by 
fitting. The dashed line and the dot-dashed line come from B = 46 uG, Ee = 4TeV, Qyeig = 0.53 


for the former and Quela = 0.735 for the later. 


be 800 Jy; while for Qyeiq = 0.735, we derive BICH* ~ 700 Jy. 


Fig. [6]shows the comparison between the electron spectra of Vela YZ with parameters in 
the previous model (solid line) and parameters in this model (dashed line for Qyeiq = 0.53 
and dot-dashed line for Qyeiq = 0.735). If we take our parameters estimated above and still 
use Vela YZ as the only predominant SNR in the fitting, excesses of e~ and e7 + e*t data 
would arise in hundreds of GeV, especially for the case Qyejqg = 0.735. Thus, we need to find 
a proper SNR as a cooperator in this case. Fig. Blindicates SNRs listed in Table [H] hardly 
have influence on electron intensity in the energy range of AMS-02, as we mentioned above. 
From Fig. Æ| we can see Monogem ring and Loop I are the potential electron contributors 
up to 1 TeV. We would like to examine Monogem ring in this model. 
ve Ld, Ring (MR) is a large shell-like structure in soft X-ray band with a diameter 

49 


~ 25° . Green does not include large X-ray regions with scales larger than 10° in his 


21 


TT rrr T1419 E rr rarr T1117 


ye 44 = 
ba —s— AMS02 TOT wr m AMS02 TOT 
re = re Monogem 
ao a — — Vela Oo MSE g 
N fo! — — MR N o — -— Secondary 
! — — Other SNRs l < E 4 
a? 85 A ETS Monogem 4 ae o 
l — l 
E> E 
O O 
q o qQ 
> => 
om oO 
oO | oO 
Yo New 
— = = 
De 9 
ui Li 
1. 
jo) 
_— 


e +e 

po T TTT T T TTT T T ft 1 2 ETEN F T T T T T J 
T = AMSO2 TOT 4 H = AMSO2 TOT 
a T VERAS =DE era 7 ro EES Monogem 
Oo eo — — Vela ee 
N Z — wr q a Me Secondary 
| — — Other SNRs 4 
ae - Monogem J E al 
1 — Secondary 
E O 
G g 
N oO 
3 

mM 
Oo | 
Yo 
— = 
5 
lu 

T 

O 

= 


E, (GeV) 


FIG. 7. Same as Fig. [5| but for Vela + Monogem Ring scenario with wela = 0.53. In the legends, 


’MR’ is the abbreviation for Monogem Ring. 


catalog hd, thus MR does not appear there. Assuming MR is in its adiabatic phase, 
Plucinsky et al. derives parameters of MR with observable quantities and its distance, 
applying Sedov-Taylor model of SNR. Plucinsky points out 300 pc should be a reasonable 
approximation for distance of MR, and consequently an age of 8.6 x 10*yr and an initial 
explosion energy 0.19 x 105! erg. Apart from this, we possess poor knowledge of MR to 
constrain its radio spectral index, total energy converted into electrons, or cut-off energy. 


We add these three quantities, &mr, Wmr and Em, to free parameters in this scenario. 


Two global fittings are done with Quela = 0.53 and Quela = 0.735, and fitting results are 
plotted in Fig. [7] and Fig. B] respectively. For the case Qwea = 0.53, the best-fit reduced 
y? is 0.398 for 180 d.o.f; if we take the latest value Qwea = 0.735, the result becomes 
x? /d.0.f = 0.396 for 180 d.o.f. Best-fit parameters are shown in Table for both cases. 
It can be seen, good fitting results to AMS-02 data are achieved when MR joins in the 
model. In e~ and e~ + et spectra, the deviation between these two results with different 


Qyela becomes clear in TeV range. 
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FIG. 8. Same as Fig. [5] but for Vela + Monogem Ring scenario with Qyeiq = 0.735. In the legends, 


’MR’ is the abbreviation for Monogem Ring. 


C. Vela YZ + Loop I (NPS) 


In this model, we largely repeat the work of the previous model, but replace MR by Loop 
I. It is well known that there are several Galactic giant radio loops potentially associating 
with SN events . Loop I, which was discovered by Large et al. and named by 
Berkhuijsen 


, is the most prominent one among these loops. Although Loop I, with a 
radius of 58°, is also not included in the Green’s catalog, it has gone through more careful 
study than MR. The center of Loop I is close to Scorpio-Centaurus OB association where SN 
events happen. North Polar Spur (NPS) is the most prominent part of Loop I, both in radio 
and X-ray maps . However, the visible X-ray contradicts with the age of 10° yr derived 
by the low expanding velocity of H I surrounding Loop I lad lag. Then a reconciliation 
is raised: Loop I is indeed an old structure, but it has been reheated by younger SN events 

, L14. The shock wave from the most recent SN event in Sco-Cen association happened 


in 2 x 10° years ago and gave rise to the X-ray feature of NPS . The comparison between 


23 


ry m ™— 7 E r T E E T T a] 
~ m 
7 —=— AMSO2 TOT a. C e~ AMSO2 TOT 
— — BKG 

a ee — — Vela (a Monsen 
N — — Loop! N —-— econdary 
Ie, — — Other SNRs 3 I QF 4 
y TERE Monogem J Y (e) Ec 
§ 3 5 
q o qQ 
> => 
om oO 
oO | oO 
Yo wa 
ao T Se 
uu Lu 

ai 

jo) 

+ 
e +e 

peo T tort T tTI T ™T 5 Ẹ T f T T T J 
ia = AMSO2 TOT J C e AMS02 TOT 
Te p= HESS = BAG ] be Monogem 
a T me VERITAS — — Vela | —-— Secondary 
N o — — Loop! E| — E 
| — — Other SNRs 4 E 
ae - Monogem J E | 
1 a Secondary r 
£ 5 
O 
N g 
3 

M 
oO | 
Yo 
Go. S 
Li 

F 

e] * | 

10 100 1000 10 100 1000 
E, (GeV) E, (GeV) 


FIG. 9. Same as Fig. [5] but for Vela + Loop I (NPS) scenario with Qyeig = 0.53. 


the continuous shell-like radio structure and the interrupted X-ray structure of Loop I also 
favors this interpretation. Thus 2 x 10° yr should be take as the age of Loop I (NPS) in 
our work, and adopting the corresponding distance to NPS—100 pc —as the distance of 
this source is more reasonable than taking that of Sco-Cen association (170 pc) as in some 
earlier researches. Similar to the previous model, free parameters associating with Loop I 


(NPS) in the following fittings are @toop, Wioop and Ee,toop: 


The fitting results are shown in Fig. [9] and Fig. for the case of Quela = 0.53 and 
Quela = 0.735, respectively. The best-fit reduced x? for the former case is 0.401 for 180 
d.o.f, while for the later case, y?/d.o.f = 0.400 for 180 d.o.f. See Table [III] for corresponding 
best-fit parameters. Comparing with the reduced x? of the previous model, the fitting effect 
of this Vela YZ + Loop I scenario has little difference with that of Vela YZ + Monogem 


Ring scenario. 
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FIG. 10. Same as Fig. B| but for Vela + Loop I (NPS) scenario with Qyeig = 0.735. 


D. Discussion 


In this section, we have proposed several source models to fit the leptonic spectra mea- 
sured by AMS-02. It can be seen that the radio spectral index of Vela YZ is crucial to the 
situation. Radio measurements from Wilkinson Microwave Anisotropy Probe (WMAP) and 
Planck may give further constraints on Qyeiq, but processed data for Vela is not available 
at present. If the radio spectral index of Vela YZ is ~0.5, other local SNRs, may not be 
necessary to join in fitting to the AMS-02 data, as in the Vela YZ model. However, if we 
choose a value of 0.735, Vela loses its predominance of the local contribution. Since the 
electron intensity contributed by the PWN is confined by e* spectrum, other local SNRs, 
such as MR and Loop I, are needed to reproduce the observed e~ spectrum. 

In the Vela YZ + Monogem Ring model, the best fitting lepton energy is 2.13 x 10“ erg, 
near the typical value of 1048 erg. However, the low explosion energy of MR derived by 
Plucinsky favors a small input energy of electrons, or a very large conversion efficiency 


of ~ 107° is required. For the Vela YZ + Loop I model, we get a larger input energy for 
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electrons of 5.08 x 1048 erg. However, this large value seems not excluded by the initial SN 
energy of Loop I of 10°? erg, which may be produced by several SN events mt Our fitting 
result Of Qjoop is 0.417, which is smaller than the value of 0.5 given by the radio observation 
to NPS between 22 and 408 MHz or between 45 and 408 MHz had. The measurements 
between 408 and 1420 MHz get a even larger value of 1.1 [117]. It is also possible that both 


MR and Loop I play important roles in the e~ spectrum. This setting may relax the ranges 


of Wmr and Aigo) Somewhat. 


IV. ELECTRON SPECTRUM ABOVE TEV 


Now we turn attention to higher energy part beyond the scope of AMS-02. DAMPE 
is expected to detect electrons/positrons in the range of 5 GeV to 10 TeV ig. Models 
proposed in the previous section have already given e~ +e* spectra extending to TeV range, 
which may be measured in future DAMPE data. However, although these models give 
different predictions in TeV range depending on the characteristic of Vela YZ, a common 
decreasing is shown in e~ +e" spectrum up to 10 TeV in each model. This means that these 


models predict no protruding structure in the highest energy range of DAMPE. 


Ground-based Cherenkov telescopes, such as HESS hal , MAGIC and VERITAS 
, have extended the measurement of e~ + et spectrum to several TeV. These measure- 
ments show that a spectral steepening appears above 1 TeV. Furthermore, HESS and the 
preliminary result of VERITAS have measured a coincident ascending of e~ + e* intensity 
in ~ 5 TeV, which may imply a feature from local sources. Unfortunately, this tendency 
comes from only the most energetic data point in both case and no measurement has been 


taken above 5 TeV. Thus it is the show time for DAMPE to examine this tendency. 


In this section, we give additional predictions to electron spectrum above TeV. We would 
study whether local sources can produce distinctive features in the high energy range covered 
by DAMPE. Although we do not intend to fit the data of HESS, MAGIC, or VERITAS 
quantitively, we need to keep in mind that the spectral steepening just beyond 1 TeV could 


be reproduced in our models. 
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A. Vela X 


There is still an important source which have not been included in our models so far: 
Vela X, as the sibling of Vela YZ. Vela X is one of the most well studied PWN powered by 
PSR B0833-45 . It has been covered by observations from radio bands to very high 
energy y-ray bands sil Lid oa} In morphology, Vela X consists of an extended halo and 
a small collimated fearture embeded in the halo, e.g., the ‘cocoon’. TeV y-ray emssion has 
been detected in the cocoon region by HESS, while Fermi-LAT observations has reported 
the presence of sub-GeV-peak y-ray emssion extending in the halo. de Jager et al. [124 
proposes a model with two populations of electrons to interpret this phenomenon: a high 
energy component concentrating on the cocoon responsible for X-ray and TeV y-ray emssion, 
and a lower energy component extending in the halo responsible for radio and GeV y-ra 
emssion. Following this idea, assuming the leptonic origin of y-ray emssion, Abdo et al. 
give a multi-wavelength fitting to the SED of Vela-X. In their results, the cut-off energy of 
the halo is only 100 GeV which is too low to help Vela-X appearing in TeV range; for the 
cocoon, although its cut-off energy is high enough, a total lepton energy of 1.5 x 10*° erg is 


too weak to produce significant structure in TeV spectrum. 


However, as mentioned above, these electron features derived_by photon emission may 
not describe those electrons we have observed. Hinton et al. provide an alternative 
model in which a serious particle escape has happened in the halo while particles in the 
cocoon are well confined. After staying in confinement status for a long time, PWN begins 
to interact with the coming reverse shock of SNR. It is the time that PWN crushes and burst- 
like injection happens. From then on, the halo has been suffering from energy-dependent 
escape, thus the cut-off energy derived by Fermi-LAT observations is such a small value. 
From some moment after the PWN crushes, the pulsar starts to inject particles to a new 
PWN, that is, the cocoon. This explains the dimness of the cocoon for its short time of 
particle injection. 

Hinton et al. gives an estimation of electron injection spectrum: E!°exp(E/6 TeV) 
with an total energy 6.8 x 1018 erg. Thus Vela X seems to produce a TeV feature. However, 
the diffusion coefficient adopted by them (1.07 x 10?7(E/1 GeV)°° cm?s~ for energy much 
larger than 1 GeV) is almost an order of magnitude smaller than ours. If we take our D(E), 


the intensity from Vela X is much larger, especially at lower energy. The spectral steepening 
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FIG. 11. Predictions to e~ + e* spectrum above TeV combining models in Sec. D] and Vela X as 
predominant contributor in TeV range. Top left and top right: Vela X combining with Vela YZ 
model; middle left and middle right: Vela X combining with Vela YZ + Monogem Ring model 
(Qvela = 0.53); bottom left and bottom right: Vela X combining with Vela + Monogem Ring model 
(Qvela = 0.735). The legends in Fig. B| Mand B]are still valid in this figure to describe corresponding 
components. Left panels show effects of varying the injection age of Vela X and right panels show 


effects of changing Do. 


indicated by HESS and VERITAS cannot be reproduced, and even the energy range of 
AMS-02 would be affected. Since the diffusion scale is proportional to \/D(£E)t, a larger 


D(E) means a faster propagation and more electrons with lower energy are able to arrive 
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at the Earth. If we want to keep the steepening feature in 1 TeV, an unreasonable young 
age of Vela X is needed. 

We combine Vela X and each model in the previous section to give new predictions. As 
discussed above, we have to keep the D(E) given by Ref. only for Vela X. Fig. 
shows the predicted e~ + e* spectra. Since there is little difference between the spectral 
shape of Vela YZ + MR model and Vela YZ + Loop I model, we only draw the former case 
as representative. The effects of varying injection time of Vela X and diffusion coefficient 
are also shown in Fig. [I] We take three different injection age of Vela X as 2 kyr, 3 kyr, 
and 4 kyr. Note these ages are still observed ages, keeping consistent with those listed in 
Table [I| With injection age increasing, the steepening structure in 1 TeV is gradually filled 
up. This requires an injection age less than 5 kyr, according with the theory that reverse 
shock has recently crushed the PWN had. Moreover, Fig. [indicates the spectrum of Vela 
X depends sensitively on diffusion coefficient (we double the Do for comparison). With the 
data of DAMPE and a clearer picture of the particle escape of PWN in the future, Vela X 


may become a potential tool to constrain diffusion coefficient in high energy. 


B. Vela YZ 


Naturally, following Vela X, the next question is if Vela YZ can play an important role 
in the highest energy part of DAMPE. The work of Ref. shows the possibility of this 
scenario after considering the release time of electrons in SNR, although they do not make a 
distinction between Vela X and Vela YZ. Erlykin and Wolfendale believe that electrons 
start their propagation after the expansion phase of SNR which hold a typical time scale of 
200 yr. This time delay is too small to change the electron spectrum of Vela YZ. Alterna- 
tively, Dorfi point out particles begin to escape from the shock front when the SNR 
dissolves in the ISM, i.e., when the velocity of the shock has dropped to the mean Alfvén 
velocity of ISM. However, even we take a very low number density of 0.01 cm~?, the derived 
release time is several times of 10* yr which is larger than the observed age of Vela YZ or the 
typical time scale of Sedov phase. To insure electrons have already escaped from Vela YZ, 
we should control its release time not longer than 10* yr, corresponding to a least observed 
age 1 kyr. The total input energy of Vela YZ is much smaller than that of Vela X, it should 
047 


be at the magnitude of 10°‘ erg. When the electron acceleration is synchrotron loss limited, 
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FIG. 12. Predictions to e~ +e? spectrum above TeV replacing the Vela YZ in Vela YZ + Monogem 
Ring model (Qyeia = 0.735) with the Vela YZ estimated in Sec. Top left: effect of varying 
observed age of Vela YZ; top right: effect of varying input energy of leptons of Vela YZ; bottom left: 
effect of varying cut-off energy of Vela YZ; bottom right: effect of varying electron spectral index 
of injection spectrum of Vela YZ. Black lines in each graph stand for total intensity corresponding 
to the Vela YZ with same line style, and solid green line in each graph represents Vela YZ adopting 
our fiducial parameters: t = lkyr, W = 2 x 10* erg, Ee = 8TeV and y = 2.0. The legends in 


Fig. Blare still valid in this figure to describe corresponding components. 


the cut-off energy decreases with the age growth of SNR. Since E. = 4 TeV is derived by 
observations of electromagnetic radiation, the injection spectrum of electrons should have 
a larger cut-off energy. Similarly, higher energy electrons should suffer heavier energy loss, 
which causes softening of spectrum with the growth of trapped time of electrons. Thus the 
spectral index of injection spectrum may smaller than that derived by radio observations. 
We set t = lkyr, W = 2 x 10“ erg, E. = 8 TeV and y = 2.0 as fiducial values and use 
some other values of these parameters to show spectral variation in Fig. In this figure, 


we apply Vela YZ + Monogem Ring model (Qyeiq = 0.735) in Sec. In this model Vela 
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YZ can bring ascending or relatively flat spectral features just below 10 TeV, as shown in 
Fig. We can imply from this figure that if the radio spectral index is taken 0.735 here, 


Vela YZ will have no chance to produce spectral feature in ~ TeV. 


C. Vela Jr. 


Besides Vela XYZ, other candidates producing prominent structure in the TeV range 
should be selected from Table Fig. shows intensity of each SNR listed in Table II, 
along with the data of HESS and VERITAS. Cygnus Loop and HB9 do not appear in the 
scope of Fig. [13] because of their very low cut-off energy. To affect the e~ + e* spectrum in 
~ TeV, Vela Jr. seems to need least parametric adjustment, namely, a larger input energy 
to leptons. Thus we examine Vela Jr. as an example. 

Vela Jr. locates in the southeastern corner of Vela on the sky map, but at a farther 
distance of 750 pc. Vela Jr. is one of those sources whose parameters are estimated _b 
fitting the multi-wavelength emission in Sec. Tanaka et al. and Lee et al. 
also conducted broadband analysis to Vela Jr. in recent years. As there have been discussions 
about magnetic field of Vela Jr., Qo and E, can be estimated by Eq. (12) and Eq. (14). We 
put aside our fitting result to parameters of Vela Jr. temporarily. Chandra has detected 
spindly filamentary structure in Vela Jr. had. This thin structure is interpreted by efficient 
synchrotron cooling of CR electrons in a strong local magnetic field of ~ 100 uG [131]. 
However, if broadband emission is modeled by leptonic scenario, a magnetic field of ~ 10 uG 
is required to explain the synchrotron to inverse Compton flux ratio. Thus we take three 
different magnetic fields of 100 uG, 10 uG and a typical value 30 uG adopted in Ref. ha. 
Here a, and B!O™ are taken to be 0.3 and 50 Jy, as given in Table [I] to calculate Qo. We 
take the upper limit of the age of Vela Jr., i.e. an observed age of 4.3 kyr. 

The predicted e~ + et spectrum of Vela Jr. combining with Vela YZ and Vela YZ 
+ Monogem Ring models (of course original Vela Jr. subtracted) are shown in Fig. 
The input leptonic energy corresponding to magnetic field of 30 uG, 100 uG and 10 uG are 
1.72 x 10^ erg, 2.80 x 10% erg and 8.92 x 10* erg, respectively. For the case of B = 100 uG, 
Vela Jr. cannot produce prominent feature in ~ TeV; when B = 10 uG, the spectral feature 
given by Vela Jr. is remarkable enough, but it needs a total leptonic energy much higher 


than the typical value of 108 erg. The result given by B = 30 uG may be what we’d like 
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FIG. 13. Electron intensities of individual local SNRs listed in Table [I| Cygnus loop and HB9 do 
not appear in the scope of this figure because of their very low cut-off energy, as shown in Table 


Data from HESS and VERITAS measurements are also shown in this figure. 


to see. However, we should note the a, of Vela Jr. is derived by the radio flux of only two 
wave bands, which is quite unreliable, and the value of a, will give a significant influence 
to the total leptonic energy of Vela Jr. with the method used here. In fact, the parameters 
given by the leptonic model of Ref. is similar to the result of our broadband fitting (the 
differences may be due to the consideration of the energy loss between injection electrons 
and radiation electrons in the work of [79]). Lee et al. also found the leptonic model 
is clearly superior to hadronic model. Thus broadband fitting may give better constraint to 


parameters of Vela Jr., which disfavor Vela Jr. as a prominent contributor in ~ TeV. 


D. Cygnus Loop 


Cygnus Loop is one of the closest SNRs to the Earth (540 pc) and is considered as an 


important local CR accelerator second only to Vela [13]. However, its steep GeV spectrum 
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FIG. 14. Predictions to e~ + e* spectrum above TeV combining models in Sec. D] and Vela Jr. 
as predominant contributor in TeV range. Top left: Vela Jr. combining with Vela YZ model; top 
right: Vela Jr. combining with Vela YZ + Monogem Ring model (Qyetq = 0.53); bottom: Vela 
Jr. combining with Vela + Monogem Ring model (avela = 0.735). The legends in Fig. [5] [7] and 
[8] are still valid in this figure to describe corresponding components. Effects of adopting different 
magnetic field are shown in these graphs. Black lines in each graph stand for the total intensity 


corresponding to the Vela Jr. with the same line style. 


observed by Fermi-LAT (see Fig. B) indicates a very low cut-off energy of electron spectrum, 
thus it cannot play a part in the models above. Cut-off energy derived by Eq. (I4) is at 
least several TeV for Cygnus Loop, so it is possible that Cygnus Loop has undergone serious 
particle escape similar to the case of Vela X. Unfortunately, although Cygnus Loop has been 
carefully studied in X-ray band with the help of XMM-Newton and Suzaku Cad kad, its 
global X-ray spectrum is not available. Thus we cannot give further inference from the 
broadband spectrum of Cygnus Loop as Vela X. Here we assume cut-off energy of TeV for 
Cygnus Loop, which is the precondition for Cygnus Loop to contribute significantly in TeV 
range. Besides, to preserve the spectral steepening in 1 TeV, considerable release time is 


needed for Cygnus Loop. We keep others parameters given by broadband fitting unchanged. 
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FIG. 15. Same as Fig. but replace Vela X with Cygnus Loop. Left panels show effects of 
varying the cut-off energy of Cygnus Loop and right panels show effects of changing the injection 


time. 


Like the former cases, we combine Cygnus Loop with models given in Sec. [HI] and also show 
the effects of varying cut-off energy and observed age of Cygnus Loop in Fig. In the left 
panels, observed age is fixed at 5 kyr, while in the right panels, cut-off energy of Cygnus 
Loop is set to be 8 TeV. 
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E. Discussion 


As discussed in the last section, it is hard to tell which local SNR dominates the electron 
excess in the energy range of AMS-02, since the spectra of local SNRs are mixed with that 
of the SNR background below 1 TeV. However, things are different above the TeV scale that 
can be covered by DAMPE. The contributions from the background SNRs may be very small 
in this energy range. If we get a distinctive spectrum above several TeV, we may determine 
the origin of these cosmic ray electrons. Of course, it is also possible that more than one 
sources shape the high energy electron spectrum, then the situation becomes complicated. 

In order to produce clear features at the spectrum beyond TeV, a small particle propaga- 
tion time and a high energy cut-off of injection electrons are essential. We have considered 
the Vela X, Vela YZ, Vela Jr. and Cygnus Loop as possible candidates contributing to cos- 
mic electrons above TeV in this section. Except for very young sources like Vela Jr., other 
three candidates need an additional release time. We find that Vela X and Vela Jr. may 
provide a very sharp rise in the spectrum at a few TeV due to the very young components, 
while the Vela YZ and Cygnus Loop may show much smoother feature at the similar energy. 
Thus with the basis of a ~ kyr injection age and a ~ TeV cut-off energy, the total leptonic 
injection energy of a source is the key factor to produce a sharp spectral structure in the 


high energy range covered by DAMPE. 


V. SUMMARY 


In this paper, we have given predictions to CR electron (plus positron) spectrum above 
~ TeV basing on the elaborate analysis of the local astrophysical sources, especially the 
local SNRs. In order to obtain a complete picture, we ensure the consistency between the 
predicted spectra and present experimental results below TeV by performing global fittings 
to all the latest leptonic AMS-02 data. We find that Vela YZ could act as the dominator just 
below TeV because of its proper age and distance. However, it should be emphasized that the 
determination of the injection spectral index of Vela YZ, which still has some uncertainties, 
is crucial to fittings. We discuss different scenarios to fit AMS-02 data corresponding to 
different values of the spectral index of Vela YZ. Other SNRs, such as Monogem Ring and 


Loop I, are also introduced, if Vela YZ dose not provide the dominant contribution to the 
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AMS-02 electron results. 

Basing on the fitting results, we discuss the parameters of several local sources, and give 
further expectations of the electron spectrum above TeV. Adopting different possible values 
for those parameters, we predict either sharp or relatively flat electron spectral features, 
comparing with the monotonic decreasing spectra of the models in Sec. All these models 
are ready for the examination at high energy CR electron detectors, such as DAMPE, which 
can reach the energy as high as ~ 10 TeV. The spectrum measurement of DAMPE and 
future anisotropy measurements |134| may reveal the origin of the high energy CR electrons. 
These results will also be important for probing the mechanism of CR acceleration in the 


sources. 
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